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Classically, black Holes have the rigid event horizon. However, quantum mechanically, the event 
horizon of black holes becomes fuzzy due to quantum fluctuations. We study Hawking radiation of 
a real scalar field from a fluctuating black hole. To quantize metric perturbations, we derive the 
quadratic action for those in the black hole background. Then, we calculate the cubic interaction 
terms in the action for the scalar field. Using these results, we obtain the spectrum of Hawking 
radiation in the presence of interaction between the scalar field and the metric. It turns out that 
the spectrum deviates from the Planck spectrum due to quantum fluctuations of the metric. 

■q" . PACS numbers: 04.70.Dy,04.62.+v 

^ ■ I. INTRODUCTION 

>>■... 

C^ ' Black holes in general relativity have shown us many fertile aspects of spacetime such as the event horizon, the 

~^ , singularity. Hawking radiation, and the connection with thermodynamics. Among these phenomena. Hawking radia- 

■ tion is particularly interesting because it gives a strong support to the thermodynamical interpretation of black holes. 

C^ \ Moreover, it reveals an aspect of quantum gravity. Usually, Hawking radiation is explained by quantizing a free field 

in a fixed black hole background [1[. To be specific, we consider a real scalar field tp in the Schwarzschild black hole 

O ' background 

^ : ds2 ^ _ (l _ !2L) dt" +{1-1^^ ' dr^ + r^ {d9'- + sin O'-dc^^) , 

where rn is a constant of integration. When the radius coordinate r becomes rn, the spacetime seems to be singular. 
Of course, we know it is just a coordinate singularity and the position r — rn is the event horizon of the spacetime. 
P*" [ The action for gravity with a scalar field reads 

QQ 1 r 1 r 
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l/~j ■ where G, g and R are the Newton constant, the determinant of the metric g^i^ and the scalar curvature. Classically, 
^^ ' nothing can come out from the inside of the horizon. Quantum mechanically, however, it has been shown that black 
^^ . holes can emit radiation, the so-called Hawking radiation [l|, y] . Remarkably, it turns out that Hawking radiation 
has the Planck spectrum. This notable feature conforms to the black hole thermodynamics. However, taking look 
at the above action, we notice that the metric should be also quantized. Then, the event horizon becomes fuzzy, 
that is, black holes are fluctuating. The purpose of this paper is to study Hawking radiation from fluctuating black 
holes. Technically, we consider the interaction between the scalar field and the metric fluctuations. In particular, it 
is interesting to see if the Planck spectrum is modified by the interaction [3|-|5| whatever small it is. Apparently, this 
effect becomes significant if the horizon radius of black holes rh is close to the Planck length £p. 

Recently, various Planck scale black holes and Hawking radiation from them have been considered 043] ■ These 
black holes can be created at the LHC or in the early universe. Indeed, black holes would be created at the LHC 
if the dimension of space-time is more than six and braneworld picture is correct [lO| . This is because the higher 
dimensional Planck mass becomes ~ ITeV which can be reached at the LHC. In other words, the mass of a black hole 
created at the LHC is of the order of the higher dimensional Planck mass. On the other hand, density fluctuations 
in the early universe could produce black holes with any mass. Black holes with lO^^g created in the early universe 
lose their masses by Hawking radiation and become Planck scale just now. The spectrum of Hawking radiation from 
these black holes might be distorted due to the nonlinear interaction between quantum fields and the metric. 

Here, we stress that, even for black holes much larger than the Planck length, the fluctuations of black holes 
might be relevant to Hawking radiation. The point is that the Planck spectrum of Hawking radiation stems from 
the exponential elongation of short wavelength quantum fluctuations which can exceed the Planck scale. Hence, the 
interaction between the scalar field and quantum mechanical fluctuations of the metric would be relevant even for the 
large black holes. 

In this paper, we calculate the effect of nonlinear interaction on Hawking radiation of a real scalar fleld from a 
4-dimensional Schwarzschild black hole. In order to achieve this aim, the canonically normalized quadratic action 
for metric perturbations and the interaction Hamiltonian are needed. Therefore, we first construct these actions by 



perturbing the Einstcin-Hilbert action with a real scalar field. After that, we give a formalism to study Hawking 
radiation from fluctuating black holes. 

The organization of this paper is as follows. In section|lll we construct the quadratic action for metric perturbations 
by perturbing Einstcin-Hilbert action in 4-diniensional Schwarzschild background. In section iHll we quantize the field 
and derive Hawking radiation. In section IIVI we present a method for treating interaction in spherical symmetric 
space-time and then calculate the cubic interaction Hamiltonian from the action of a real scalar field. Using these 
results, we calculate the effect of nonlinear interaction on Hawking radiation of a real scalar field. Section|V]is devoted 
to the conclusion. Technical details can be found in Appendixes. 

In this paper, we use the convention IGttG ~ 1 and write the Schwarzschild metric as 

ds^ = -f{r)dt^ + r\r)dr^ + r^jabdx''dx'' , f{r) = 1 - — , 

where jab represents the metric of the sphere S^ and M corresponds to the mass of the black hole. 

II. METRIC PERTURBATIONS IN BLACK HOLE BACKGROUND 

In this section, we derive the quadratic action for metric perturbations h^i, in the Schwarzschild black hole back- 
ground gliJ . Note that our derivation is slightly different from the original one by Moncrief [Il|. Now, the metric is 
given by 
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Here, we can use the general coordinate invariance to reduce the number of variables. Under infinitesimal coordinate 



transformations 



x^ + ^^ , hfj^i, transforms as 



'^/ii/ ^ ^^.u I s^'M ' s?^;/^ 5 
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which shows that there arc 4 gauge degrees of freedom. Hence, we can fix 4 variables. Moreover, we can classify the 
remaining 6 variables into 4 scalar type and 2 vector type ones. Once the gauge is completely fixed, we can obtain 
the quadratic action. After eliminating unphysical variables using constraint equations, we finally obtain the reduced 
action for physical variables. 

It is easy to extend the analysis in this section to more general cases. In Appendix |B1 we present the analysis for 
topological Schwarzschild black holes with a cosmological constant. 



A. Scalar Perturbation 



First, we consider the scalar type perturbations 



/ fH H, 
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where \a denotes a covariant derivative with respect to jab- Under the scalar type gauge transformations ^^ with 

Cm= (^0 ,Cr ,a|a) , (4) 

the perturbed variables change as 

H^H + 2io/f-f'^r, H,^H,+S.'o+i-.f'^o/.f . v ^ v + Co + ^l , (5) 

H^H + 2fi + f'Cr, iv^w + Cr+£.'L-^^L/r , K^K + 2fCr/r, B^B + 2Cl, (6) 

where ' and • denote derivatives with respect to r and f, respectively. Then, choosing ^^ = —B/2, ^o = —v — B /2 
and ^r = ~^A'/(2/), we get the following metric perturbations 
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Using the above gauge, we obtain the quadratic action 



e,m>o 



drdt 



^\Hum\ - 7sW*g„iHiim + y\'^^^\ + 2r-i?^*,„-ffif 



^ / drdt 



V + r.f , 



ilr) 



2^ -wlmHhnA :^\Hhn\ + 7s ^ k«m | ~ ^Hl^q2hn 



{C.C) 



HuQ^ - jsWioHuo + yw^o + ^rHioHn 



V + rf 

-Is — W£0 Hi, 



ifr 



f 



1 



2r "" '" 2 

where * means complex conjugate and we have defined 

qiim = 2rfHer,i - 2^sfwi>m , <72£m = 7s f Him + / Wl 

It is easy to solve as 



Hio + 7s^w^o - o^ffi92ffl 
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where we have defined T{r) = rf — 2/ + 7s. In Eqs.®, ©, and pl?|. we expanded metric perturbations in terms of 
spherical harmonics Yim like as -ff = ^ HirnYim- Here, 7^ = (.{(. + 1). Since spherical harmonics satisfy the relation 
^£m = (~l)™^£'^m' metric perturbations Hi„i satisfy reality conditions _ff^,„ = (~l)'"-f^i'-m- This allows us to relate 
negative m modes to positive m modes. Hence, we can consider only tti > modes in ([5|). Note that m = Q mode is 
real but to > modes are complex. 

From now on, we concentrate on ?Ti = mode and omit the indices. We can get the same result for to > 
modes with similar calculations. Taking the variation of the quadratic action ([8]) with respect to ifi, we obtain the 
momentum constraint equation 



Hi = w 
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This equation enables us to delete Hi from the action ([8]). Substituting this result ([TT]) into the action (j8]), we obtain 
the action 



drdt 
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H^ + 2rHw - 7s ^' ^^ ' wH + ^H^+j^J^w^ - -Hq2 
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Taking the variation of the action (|12p with respect to the Lagrange multiplier H, we can derive the Hamiltonian 
constraint 



92 = 



(13) 



Then, substituting q2 = into the action (|12p and using the relations ([T0| . we can express the action by qi. With a 
new variable tp^ defined by 



r 

the quadratic action takes the canonical form 
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where r* represents a tortoise coordinate and we have defined the potential function as 



Vz = 
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4(r2/' -2rf + 2rX + 2r) 



8A(1 + A)^ + 4(1 + Xyrf + 2(2 + A)H/ ^ + r'f 



2 J-'2 I ^3 f'3 



+4/2(2A + rf') - 4/ {2A(2 + A) + 2(1 + X)rf' + r^f'^^ 
2/ A2(A + l)r3 + SA^A/r^ + 9\rAP + 9M'^ 
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We can also calculate the quadratic actions for m > modes. Thus, we have 



(16) 



J2 ^jdr*dt {%^U^ - \dr^i^L? - Vz{rMM } + \jdr*dt [{dt^f,f - {dr^^^f - ^z(r)Vfo'] , (17) 



which gives the famous Zerilli equationjl/ 



5,VL + dl,^L - Vz^L = 



(18) 



For later purposes, we deduce the relation between ipf^ and metric perturbations. From Eqs.(j9l), ([T0|. ([TT|) and 
14]), it is easy to get the relation between -0^ and metric perturbations except for H. The results are 
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In order to get the relation between tpf^ and Him, we need to look at the equation of motion for metric perturbations. 
Taking the variation of P^ with respect to H and w, we can get two equations 
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Him - 2rWtra " Is 7, ^^m + (^/) ^Im - -^Him + r f H ^^ = , 

7s 2r 2 

-2rHim ~ 7s t; Him + 2-is^wim - TT (/ )Him - IsJH^.^ = 

2r r^ 2 r 



Eliminating i?£„j from the above equations, we obtain the equation for Him 
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Using this equation, the relation (J19p and Zerilli equation, we can show that Him is related to Him and wim. as 

Him = Him — {2fwim) ■ (23) 

Since Him and wim are related to ipfm by Eqs. ([T9l) - this equation implies a relation between Him and tpfm- 

B. Vector Perturbation 



Next, we consider the quadratic action for vector type perturbations 



flni/ 



sym sym 



I'a 

Wn 



2 [C^a\b + Cb\a\ 



(24) 



where Va, Wa and Ca satisfy vj°' — wj'^ = Cj'^ = 0. For vector type perturbations ([24]), the gauge transformations 
with 



read 



eM = (00,Ca) , 



Va -^ Va+L , Wa -)■ Wa + ^a - ^iaJT , C a ^ C a '\^ 2^^ 



(25) 
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where ^a satisfies the transverse condition. Then, we can completely fix the gauge by choosing ^a = —Cal'^-, which is 
called Regge- Wheeler gauge. In the Regge- Wheeler gauge, metric perturbations are written as 



Wa 

zwa 
sym sym 
sym sym . 

Using the metric perturbations (j27p . we can calculate the quadratic action 
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where metric perturbations are expanded by vector harmonics V^°^ like as w"- = J^wimV^-f^. Since vector harmonics 
satisfy relations V^,^ = (— 1)'"141*„, wim satisfy reality conditions wim = {—^)™'we-m- Hence, we can consider only 
m> modes in (f28| . 

For a while, we concentrate on tti = mode and omit the indices. Taking the variation of the action (|28p with 
respect to v, we can derive momentum constraint equation 



r (ui H — V — V ) 
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We have to solve this constraint equation. However, it is not easy to do so. In order to circumvent this difficulty, we 
move on to Hamiltonian formalism. From the quadratic action ([28]) . we can read off the Lagrangian 



L = (2A + 2) / dr 
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and also define the momentum conjugate to w as 



p= -T = {2X + 2){w + -v-v 
ow \ r 
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Then, the Hamiltonian H is given by 
H = I drpw — L 



= (2A + 2) I dr 
Now, the constraint equation 

gives the relation 
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Substituting this relation into the Hamiltonian ((32)) . we obtain the reduced Hamiltonian 
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2 2A XJ \2X + 2J 4A V2A + 2, 
We further perforin the canonical transformation [w^p] — > (Q, P) to interchange the role of w and p 

Q = -P , w^ P , 
which leads to a new Hamiltonian 
A'(0,P) = H{P,-Q) 

= (2A + 2) /"dr 
Using the Hamilton equation for Q 
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we can go back to the Lagrangian formalism 
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With a change of normalization 

the quadratic action (|57| becomes 
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where r* represents a tortoise coordinate and we have defined the potential function 
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So far we have only considered the m = mode in ([28|) . Of course, we can also calculate the action for m > 
modes. The result reads 
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which gives the Regge- Wheeler equation [l3| 
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Finally, we clarify the relation between ipf^ and metric perturbations. From Eqs. ([33| . ([35]), ([36| and ([38| . the 



relations are given by 
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III. HAWKING RADIATION 



111 this section, we review Hawking radiation using a spherically symmetric collapsing model[l|, |lj]. In this model, 
it is assumed that outside of the collapsing matter can be described by the Schwarzschild metric 



ds' = -f{r)dt' + l/f{r)dr' + r'-fabdx^'dx'' , 



(44) 



where /(?■) = 1 — 2M/r. Here, M and jab are the ADM mass and the metric of 5^, respectively. The action for the 
real massless scalar field 
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can be written as 
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where r* is a tortoise coordinate. Here, we used the harmonic expansion </? = X) 4'hnYim/T- Note that (p is real. 
Hence, complex coefficients (jy^rn satisfy reality conditions 4>irn = (^ 1)™0^_,„- Then, the equation of motion outside 
the matter (p''';p = can be written as 



a^,/),™ - a2.0,,„ + \/^0,„ = , T/^ = /(r) ( L + :^^^-tll 
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It should be noted that the quantization of metric perturbations can be performed in the same way as that of the 
scalar field explained below. 



A. Quantization 



Let us quantize the scalar field ^ by promoting Lp to the operator: 



<fi 






au,tm.-YtmU^t{t,r*) + a'^f^^-Yl^ult{t,r*) 
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In the above formula, -^=^Y(rnUu>f. are complete orthonormal basis satisfying Klein-Gordon equation 



d^^u^i + dl.u^i - fir) I L + ^M±}1 U„, = . 
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Note that this equation is independent of tti, so we used Uuji instead of Uajgm- Here, the norm is defined by the inner 
product 






(50) 



where S is a Cauchy surface and dY,'^ means area element times a unit normal vector. Now, the canonical commutation 
relation implies the following commutation relations; 



Using the property of spherical harmonics Y^m = (— l)™i^/_„j, we can rewrite (|48| as 



f^^j '^~7k= aa;fmUa;^(i,r-*) + (- l)''"aL-m<^(*' ''* ) 



Yt 



r 



(51) 



(52) 



Let us define (j)£rn by (/3 = ^ 4>imXhn/f- Then, we have 

1 



aim 



= duj- 



bujemu^eit,r*) + dli,^ulg{t,r*) . (53) 



/2uj 
From Ea. ((5^ . we see the following relations 

buieTTi = aujim, d^tm = (— l)'"aa;f-m • (54) 

These relations are equivalent to 

'/'L = {-^r^i-m. , (55) 

which is nothing but the reality condition for ip. And this reality condition suggests that we can only consider m > 
modes. From (j5ip and (j54p . we can deduce the commutation relations for b^im and d^jim for to > modes. The 
results are as follows. 

• TO = modes 

We can define fe^j^o = O'ujIo = dujmi^ '^'^i) from the relations (|54p . Using the commutation relations (|5ip. we can 
show that 

[Ca;f,C^^,] ^ S{ll!~Ll!) 5^^' . (56) 

The commutation relations show that 0^0 are real scalar fields whose annihilation operators are c^ii- Note that 
these scalar fields dm have the action (l46l). 



• ?7i > modes 
From the commutation relations (j5ip . it is easy to show 

Thus, (f>im {m > 0) are complex scalar fields whose annihilation operators are bi^im and dujtm- Note that complex 
scalar fields (pim have the action _ 

B. Bogolubov Transformations 

To introduce the Bogolubov transformation, we define complete orthonormal basis 

VuiUm = -^==-YlmUul ■ (58) 

V2ujr 



Suppose that there are two kind of complete orthonormal basis of (p, {vuigm} and {v^e„i}- As both {vi^im} and {vi^im} 
are complete basis, we can expand ip in two ways 

im £m 

where a^im and (ii^im are annihilation operators for v^im and Vi^irm respectively. Using these annihilation operators, 
we can define two kind of vacua as 

a^im\0) = , a^em\0) = . (60) 

Due to the completeness of {vi^i„i} and {wi^^m}, we can expand Vujim as 
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where Qf^^^^^'^'^' and P^im ui'i'vi' ^^^ called Bogolubov coefficients. From the above formula, we can deduce the 
Bogolubov coefficients as 
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Since both complete basis are proportional to the spherical harmonics Yem, Bogolubov coefficients are proportional 
to Sf_j_^ and S^> , that is 
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Substituting (J6T)) and (|62|1 into (|59|) and using (J64|l . we can get Bogolubov transformations between auitm and a^ 
as follows; 
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Furthermore, substituting the second relation in (pSl) into the first relation in (p5|) . we see that Bogolubov coefficients 
satisfy 
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So far, we have considered Bogolubov transformation for ip. We can also describe Bogolubov transformation in terms 
of (j)im {m > 0). In order to get Bogolubov transformation for c„£ which is defined by c^i = ai^iQ, we set m = in 
(p5)) . Then, we get 

C"^ = J doj' [a^' e,ujec^' I + l^l'i,^f\^ . (67) 

Similarly, setting ttt, > in (j65p . we have Bogolubov transformation for bcuim which is defined in (|54p : 

And we obtain Bogolubov transformation for di^im by considering tti < in (|65p 



d<^fm =- I duj' { a^' t^^td^' irn + 1^1;' e,u,A im 



Finally, we calculate particles created due to the time evolution of space-time. It is assumed that {vojem} is a 
natural base at the initial time ( "in region") and {vi^im} is natural at the final time ( "out region"). Furthermore, 
we assume that the initial state is |0). Note that this state does not change because we use a Heisenberg picture. 
Under these assumptions, the number spectrum of (p particles at the final time can be expressed by 



J2m^i,n\0) = J2^0\al,^a^em\0) 



J 



i 0<m<f 

Then, using Bogolubov transformations for h^gmi c^i and d^^im, we can calculate the number spectrum 



E(oi^-^™|5) = E 
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(71) 



This shows that we need /3^^ j ^ defined in (|63p and (|64p in order to know the number spectrum. 
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FIG. 1: Penrose diagram for a collapsing model. In this figure, A means an affine parameter 



C. Hawking Radiation 



Let us consider Hawking radiation by using spherically symmetric collapsing model described by the Penrose 
diagram in FIGHJ In this figure, /~ corresponds to "in region" and /"*" represents "out region". 

In order to calculate /3^^ j g, we must define natural basis in "in region" and "out region". We assume that outside 
of the collapsing matter is described by the Schwarzschild metric. Then, 0£m satisfy 



dti'tm - d^,(l)im = 



(72) 



in this region. Here, V^ is ignored because V^ dumps as 0(l/r^). Noticing ([72)) . we choose the natural basis {w^} 
and {u°"*} defined by 



-> e 



-iLj{t+r*) 



(near / 



-iuj{t—r*) 



(near I'^) 



(73) 



Given the complete basis, we can calculate P^g ^'^ defined in ([63|) and (|64| 
We have to choose S to calculate the inner product in the formula (|55)) 



Here, we take / . Hence, it is necessary 



to know the behavior of w°" near /~. For that purpose, we use geometric optics approximation where the worldline 
of a particle is a null ray of a constant phase and trace this null line backward to I^ . This approximation is quite 
good near the horizon H^ because the frequency effectively becomes very large in this region. Hereafter, we write 
u = t — r* and v = t -\- r* . Furthermore, we define A as an affine parameter with respect to n^ which satisfies A = 
on "jh- Near iJ^, this parameter can be expressed as 

A = — Ce"**" (near horizon) 
A = Dv (near H~) ' 

where k = 1/4M is the surface gravity and C and D are some constants. Because A — const, line is also null line, 
geometric optics approximation suggests that A = const, surface is equivalent to a constant phase surface. Then, the 
constant phase surface A = — e satisfies 



e = — Ce "" (near horizon) 
-e = Dv (near H~) , 



(74) 



and so 



. = -ilog(-§. 



(75) 
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This means that u°"* approximately behaves as 

out / (^°^ ^ > '^) 

"^'^ - j 1 e't^iosC-*^) (for w < and |w| < 1) ^''^^ 



near /~. Then, we can calculate the inner products in (|63p using (|76p . Note that ([76]) is a good approximation in 
|w| <C 1 but the effect from this region is dominant because the phase of u°'^* diverges. Therefore it is also a good 
approximation to use the expression ([75)) on /~ in calculating inner products. The results of calculations are 



^fj-c.,[^H-Dv/C))d. 



^\ " " g^-/2«(^y^^/Kp i"^ ^ ^^1 (77) 



and 



\l^£j--e.,(^^H-D./C))d. 



This gives 






«...'/ =e2-"/«|/3^,,,r (79) 



Using the relation 

d^' (|«C.£V/ - l/^c..,.'/) = 1 (80) 

derived from ([66]). we obtain Bogolubov coefficients /? 

d^'lP.i.^'f = ^.J. - 1 ■ ^^1) 

Substituting this result into ()7ip . we find that the observer in the asymptotic region see the thermal spectrum of Lp 
particles. 

Actually, Ea.([8T|) diverges and is proportional to the delta function (5(0). We interpret this as 

1 I^TVL 



(5(0)= lim — / expi(wi-w)i (82) 

and introduce the regularization (5(0) = L hereafter. Namely, what we obtained is the number of particles created per 
frequency and per unit time. In the above formula, we have discarded this factor. However, we should keep in mind 
this fact in the following calculations. 

IV. HAWKING RADIATION FROM FLUCTUATING BLACK HOLES 

In this section, we consider interaction between the scalar field (j) ^-^d gravitons ip. First, we introduce a method 
for treating interaction in the spherically symmetric space-time |14{ . Next, we calculate cubic interaction terms from 
the action ([45]) . Finally, we evaluate the effect of interaction on Hawking radiation. 

Here, we follow the notations in Sec lIIIi namely, we use (pim defined hy Lp — ^ 4'imYirn/T a-nd annihilation operators 
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A. Interacting Quantum Fields in Curved Spacetime 

Let us consider the Lagrangian 

•^ = -Co + C-mt , 

where Co represents the free field part of the aetion and Cint is the interaction part. Furthermore, we assume that Cint 
vanish both in the remote past (or "in region") and in the remote future (or "out region"). Under these assumptions, 
we can define asymptotic fields 



ip —J- ip^" (for in region) , ip — > 1^9°"* (for out region) , 



(83) 



where (/3™ and (p°"* satisfy equation of motion derived from Cq. Therefore, using basis {-^-^u™^{t,r*)} in "in 
region" and {-h=YimU'^/{t,r*)} in "out region", we can expand asymptotic fields v?™ as 



y^ 



Y, dujia, 



1 ^«m in , in t ^ ^em tn . 



2w r 
1 Y, 



2u r 



2uj r 



'tn out j^ -in t tm „,out^ 



2uj r 



(84) 



and 1^°"* as 



V> 



Yd. 

era •' 



nut ^f-rn ^in .nut f ^tm „,in' 



2u} r 



2uj r 



ut 
'ujlm 






From the above expansion, we can define four kind of vacua 

aL"™|Oin)-0, CfmlO in) - , a^^^lO out) = , a™,L|0 out) = 



(85) 



(86) 



Here, |0 in) is a natural vacuum in "in region", while |0 out) is natural in "out region". Other vacua are introduced 
merely for the computational purpose. 

Now we calculate the number of particles created due to the time evolution of space-time in the presence of 
interaction. Given the natural vacuum |0 in), the number spectrum of Lp becomes 



^(0 in|iV°^*jO in) = ^(0 in|a°';:iC"™|0 in) 



tm 



Im 



o<m<e 



= E[(0i^lC"*^^"^"*|0i^)+ E {(OinrGrCL|Oin) + (Oin|J°t!rfr.L|Oin)}] , (87) 

Tout-fTout In ;^\ T^It^ T3^n-^1,,K^,, f ,.r,„e,f^^,Tn„ + ;^„o /[7T^ f^^ Znut 



First, we evaluate (0 in|6°"Jj6°"*„|0 in). The Bogolubov transformations (p5|) for &°"^ read 



lOUt 



' ujt.uj t uj tin 



/3% ,„d°V* A . 



Then, b°Jg^b^f„^ can be expressed by 



lOUt^lOUt 

"uj£m"ujem 



dio' 



+ / duj' / duj 



7 OMtf lOUt 



out-f mut 



* youi-^ inut 1/3 Q* jouzj inu 

''oj£,uj" £"^u:£,u:' £"u" £m uj' £rn "*" Poj£,t^' £Pu:£,u:" £"'uj" £m ui' £m 



-(ot, ,f , ," fP* e '«&°"1 rf°"/ + h.c) 
\ ia)(,iaJ t:'^uj£,uj £ u £m u £m ' 



(89) 
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The expectation value of this number operator in the state |0 in) is given by 

2 



(0 in|G„lC.L|0 in) = <ko' 



P. 



+ dio' duj" 



<^u,fu,"f(^\ '.(Oiii|&°"*l &°"*, |Oin) 



~P.,p.,'pP% ,v(OinM™*! d°'^* 10 in) 

-^^la^tu"iP% '.(Oin|&°"*l d^f/lOin)) 



(90) 



The first term of the right hand side of this equation is the same as ([7T|) and so this term represents the conventional 
Hawking radiation. Hence, other terms in (|90p are induced by interaction. In fact, these terms vanish if there is no 
interaction. Then, in order to know the effect of interaction on Hawking radiation, what we have to do is to calculate 
these terms. For later purpose, we decompose these terms into u = u cases and lo ^ w cases. Using L = (5(a; = 0), 

for example, f duj' f dio" a, ,e , ," ^a* „ ,„(0 in|6°?,*J 6°"* 10 in) in (IMl) can be written as 

-^ ' -^ ^ (^t,uj t uji.Lj i^ ' UJ im UJ im' ' ^ '' 



OUtf lOUt 



duj duj a. >> ,a\ ,„(0 in|6°r! 6°? 10 in) 



= I duj du: a^e.uj"iali.uj'i 



(0in|&™7 6°?* lOin)-^^^^:^ ^(o in|6™*J 6°?* |0 in) 



L 



- '^^^'~"^"^ 0in|b°?/J 6°^* lOin) 

T ^ ' UJ im !»J im. I ' 



du 



(0 in|&™*; 6°V* |0 in) 

^ ' UJ im w im' I 



duj duj a^ij^^Q*^^ j^ 



(0 in|6T/J h°% |0 in) - ^^'^ "^ '^ (q inj^r/ fe""^* |0 in) 

^ I UJ im. w im, I ' T \ I tj im. uj im. ' ' 



(91) 



The second term in the above equation reads 



■ ^,out^ ,out in ;„\ ^(^'-^")/n =„i;,o«*t ;,o«t m ;„x _ J (0 i" I &7*L ^°/ L I « i^) (for ^' ^ ^") 



(0in|6TI 6°^ |Oin 



^ -(Oin|&°?7 6°^ |Oin) 

^ ^ ' LJ i7n UJ im' I 



(for OJ = CJ ) 



(92) 



Therefore, Eg . ([911) can be rewritten as 



- /do; 



^ujt,uj'i 



(0 in|6™*; 6°'^* |0 in) 



+ / dcu duj 



0^uji,uj" e<^uji,uj' i 



% v(Oin|fo°"f; fo°' |0 

uj£.UJ r. ^ '1,1 Prn. UJ r.m' 



UJ £rn UJ hn I 



(93) 



Similarly, we can evaluate expectation values (0 in|d°"„jd°"^|0 in) and (0 in|c°'^ c2"*|0 in) in (|57|) as 



(0 m|C™C^L|0 in) = / dJ 



(0 in|C7'c°f |0 in) = / du:' 



doj' duj" a.,fn.a\ ,,(Oin|d™;j d°V* 10 in) 

' '■■'t.uj I ujI.uj i^ I UJ £m uj tm' ' 

+Pujiuj'iP\ ",(0in|6°f*^ 6°V* |Oin) 

i^ujt,uj tr u)£^u) £\ I ^ ^j^ tj ^rn I ' 

-23fi (a . .-,/?% ,,(0 in|d°f,*J 6°ff |0 in)| 

1^ i^<;,iA; ti-'uji,uj i^ ' UJ im UJ im' 'J 



rfw' / rfw" 



(«^.,^".<...'. + /3..,.'./3:,..",)(0 in|c:"*!<'^|0 "^) 



-25ft{a.,v'./?:..'.(o-i^::^^:"'ioin)} 



(94) 
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These results are almost the same as (0 in|6°"J^&°£^|0 in) and the terms other than the first term of the above formula 
have information of interaction. Substituting ([TO)) and ([M]) into (|57)) , we find the deviation from the Planck spectrum 
due to the interaction as 



E 



duj' / duj" 



(^.i,."i<,W, + Pu.Lu.',PU^"A (0 in|c™*lc°V*|0 in) 



+ y f(0in|6°?/; 6°V* |0 in) + (0 in|d™/; d°V* |0 in) 

0<m<f 



0<m<l 



(auiu"iP% 'J y" ((Oin|6™/; d°f; |0 in) + (0 in|d°f/; 6°?*^ |0 



(95) 



In order to calculate (0 in|6°" „ 6°"* 10 in) etc., it is useful to define the basis 



|nt • • • ; ?ic • • • ; »^d • • • out) 



(&°"*t)"' . . . (c°"*t)"^ . . . (do»«t)"'' . . . |o out) , 



y/ubl ■ ■ ■ Ucl ■ ■ ■ nj ■ 
which satisfy the completeness relation 

|n& ■ ■ ■ ; ?^c ■ ■ ■ ; "d ■ ■ ■ out) {rib- ■■■,nc---;nd-- ■ out| 
(ni,- ■■;nc-- ■■,nd-- ■ out|7ib ■ ■ ■ -.n^- ■ ■ ]nd- ■ ■ out) 



1-E 



(96) 



Substituting (|96|) into just after (0 in| and just before |0 in), for example, we can calculate (0 in|6°" J &°V* |0 in) as 
follows: 

• If OJ = 10 , 



(0 -i^:'i^^'Lio in) = y: e ".^v™- 

"bu.'«m {ot/iers} 



-^ 1(0 'Mn^jt^] {others) out)| 
"6c^'«m;{o^/iers} out|nj^'(,„; {oi/iers} out) 



(97) 



where i = 5{lu = 0). 
• If oj 7^ oj , 



(0 -i^:;^„^^'Lio in) = E E E 



L 



■^ , ■'^-^ , (n' ,„ ,n'' „„ \ {others} outln' ,„ , n" „„ : {others} out) 



\/<^"to«Kc.'£m + l)(0 in|n^^,,„^,n|,\,,„^;{ot/iers} out) 



X (^fi,y^^. + 1' "(,„/>„,, - 1; {others} out|0 in) 



(98) 



CJ £ CJ 



Similar rules apply to (0 in|c°?, Jc°'/*|0 in) and (0 inlcf"" J d""* 10 in). In the same way, (0 in|c°» Jc°?„ 10 in) becomes 

• If W = UJ , 



(0in|c:?*tc-*t|0in) ^ ^ E ^^ i^c.' i + ^)i^c.'i + 2) 

"ctd'f {others} 



(0 



"-cuj e 



2; {others} oiit){n^^'f, {others} out|0 in) 



{n^^'f, {others} out\n^^' f, {others} out) 



(99) 
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If w 7^ w 



(0in|c:^c-t|0in) = E E E 



L 



, (n ,.,n" ,, -i others} outln' ,.,n" „ „:i others} out) 

n" „ {others} ^ cuj i' clj i' '- ' ' "■' "' "•' '' >- J ' 



CCJ £ cu 



("' ,V + !)("" ,". + 1)(0 in|« ,v + 1' ^ru,"e + 1; {ot/jers} out) 



cu I 



cu e 



'^('^cuj'v'^cw"e'^^°^^^'^^'^ °"*|0 ^") 



(100) 



And also (0 in|d°"*/ b""!} 10 in) is given by 

(0 inM-l6:^jO in) = E E E 



71 , n ,, 

duj £m buj £n 



r , . \^'^ 'n )^r ",; 'Aoth^Ts) outln' , ,n' ,. -Aothers) out) 

{others) ^ dtu £m' bu) im^ "- -* ' dcj £m' 6cj ^m' "- -* ' 



\/K^'^m + l)("L"£™ + l)<'^ in|ji;;^,,^ + l,7iL",„ + 1; {oi/iers} out) 



x("rfM'^m'"h<.'Vm;{o^^ers} out|0 in) 



(101) 



These results will be used in the calculation of S-niatrix. Before moving on to calculations of S-niatrix elements and 
clarify the effect of interaction on Hawking radiation, we need to calculate the interaction Hamiltonian. 



B. Cubic Interaction terms 



In this subsection, we derive cubic interaction terms in the following action 

1 



d^J^q'"' 



99 'P,f^f,u 



We put the metric 



JO) 



5/. 



- 0^°) + h 



where gp_y is the background Schwarzschild metric. In the Regge- Wheeler and the convenient gauge, /i^^ reads 



'^tiv 



/ /5 Ri 

sym H/ f 


Va \ 

Wa +W\a 


sym sym 
y sym sym 




0/ 



Using this gauge, we obtain the cubic action 



- [r'Vl\{H-H) 



-■- -2 I J- '2 , 7 



We see that scalar and vector perturbations are decoupled, so we will calculate them separately. 



(102) 
(103) 

(104) 



(105) 



Vector perturbations 



In the action (|105p . terms including vector perturbations are 



d'^Xy^ 



--V''ipip\a + fw"'^ V\a 



(106) 
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We expand these variables by harmonics; for example, 



Im 



WtmV[ra ) </'=/] ^imXlm 



im 



where V^^j are vector harmonics and Yim are spherical harmonics. The result is 



D", 



£m 



m :t m 



drdt 



-j'>^em'Pi',n"Pi"m" + f ^Inif f m' '^ i" m" 



where 



£ rrt ■£ m- 



dQd^^VZY,,„.Y, 



(. m \a 



(107) 



(108) 



It is easy to show the relation D^™ , ,, ,, = — D^,™ „ , , using the property Vf°j|a — 0. We also sec that 
Dp" , ,, „ are pure imaginary numbers proportional to <5„ „'+,„" because both vector harmonics and spherical 
harmonics are proportional to e*™*. Furthermore, because DJV' , ,, ,, are pure imaginary and both harmon- 
ics have property such as Yim = (— 1)"'!^"^, Dp^ , „ „, we can derive the following formula. Dp" , „ „ = 

= (-l)"+"'+'""+'^/"_™'^,"_^"- Note that we omitted EemEe' m' Ei" m" in (HHB. Thus, using 
- (f'im/f, we can get cubic interaction terms 



H rrt \t m 

(l43t and tpf„ 



DT , ,n n \ dr*dt 






+ -<r*'/'£'™'5.*' 



(109) 



2. Scalar perturbations 

Terms including scalar perturbations in the action (jlOSp read 

1 



1 /■ o ^1 



2/ ^'^^^H-H) 



2 '2 I 

-ip + fifi +—^(pia'P\b 



1 



r ^/7 



H 

1 



■tp^ + 2Hnp (p - fH(f 



'2 



2f I 



w'"V3 ip\a 



Using expansion by spherical harmonics like 

we can write the action in the following form 
1 



<^ 



= 22 ^t-n^^t 



e.m 



o i m -.i m 



drdtr 



J [ifWinS - Hl„}j (fi'^'iPg"^" + f {[fwlJ - ij;,„j 'P^'^''f\ 



m ' I m 



2 w 



2rH* 



Im 



\a 



■'^fKmfe'm''^tm\a 



where we have defined 



(110) 



(111) 



Cpl,."n." ^ I V^YLy,„.'Ytm" ■ (112) 



These coefficients CJ^ , ^,, ,, have almost the same property as D^J" ' «" " ■ It is easy to derive the relation 
'^/" './'™" = C'/'L"./™'- Because of Y«„ ex e™*, Cp"^, ,^„ ^, are real numbers proportional to ^™,„'+„". The 

^™, „ = Cf,™, ,„ „* = (_i)™+™'+m"cf-™, ,,_ Using (HSl), ^ and 



m :t m 



reality of C!/" , „„ „ implies that CS' 
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ftm — 4'im/f, wc can derive cubic interaction terms for scalar perturbations 

2r 






___a2.V>f„: + ( ^ + Mr) j dr.^l: + B{r)i;f,: J ^,,^4," 



IsJ V 7s/ / 



/2 



2/ 



'?^£'m'^^ 



m'-r-V^/V" +dr'(l)e'^'9r'<l>i"m'' 



where 



+ / ^a„.d£_ / ^ - E(r) J ,1,^A 0,.,„, (^-^#,.,„. +3,.*,. 



7s 7s-r 7s V 7s-r / 7s 7s i 

, . T' 2 rfr' 2rf (t'\ ,, 4(7s - 2) 
D(r) = — + -{ — ] , E(r) = -^ . 



(113) 



(114) 



These function A{r) etc. are regular outside the horizon. However, the other coupling functions in (|113p are singular 
at the horizon because these terms are proportional to l//(r). We also have a relation 



A{r)+C{r)+E{r) =0 . 



(115) 



Note that we ignored the terms which is related to equation of motion for ip in (|113p . This is because these terms do 
not affect our calculation. 



3. Total cubic interaction 



We notice that the cubic interaction Hamiltonian have the following structure 

H^nt = C*/™';/'™" Yl / dr* dthM{r)ipf^ 

irit 

y^ / dr*dtwint{r)i^. 



m Tii m 



~^-^e'm';e"m" 



RW* I 7 

im Ve'm Ve"m 



From the interaction Hamiltonian pi6p . the first term of (|109p can be reproduced by the replacement 

If • - 

v/A(A + l)r^ 



(116) 



(117) 



Note that m takes the integer in the range —i < m < £; m and m take the value in the same range. As is shown in 
Sec lIIIi (f)im (m > 0) are complex scalar fields and (fiiQ are real scalar field. By using reality conditions, we can rewrite 
the interaction Hamiltonian (|116l) using only m > modes. From now on, m moves in the range < m < i unless 
we specify its range. With this notation, each term for scalar perturbation in (J116p can be rewritten as 

A^ A^ A^ Z^ ^l'm';l"m"'^in 

1,1' , I" ^<™<^ I' <m' <l' I" <m" <l" 



+(-1)™ ci 



£ rn 

A' _i_ il^.^A . . JJ 



• i^iO(t)e'm' 4'1"'' + ^m<Pt'm' V„," ) 



+(-1)™ c 



£m 



{^ern4'l,(l)i"rr," + V'fm<?!>L , </>/' m" ) 



(118) 
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graviton 
(©",1") 




(00,1) 



^((oM') 



FIG. 2: This diagram afTects S-matrix we consider. This diagram corresponds to K^ 



graviton 
(COM') 




^ (CO, 1) 



FIG. 3: Loop diagram which appear in the first order calculation. 



z ^^A A.. ^^ Note that there 
However, such terms must vanish due to the faet 



In this deeomposition, we use the property CP,.,, ,, and reality conditions for V'lm and 
exist other terms in ^^ such as C^^/°„,.^„^„'(/'fo'^«'m"^£"m 

f.' m' -,1" m" "m,m +m ■ 

Similarly, the terms related to vector perturbations read 



£ m i m 






-(-1)" D' 



m :e —m 



■{ipi 



V'm" -V'«m0t .0/V") 



(119) 



C. The spectrum of Hawking Radiation 



In this subsection, we calculate the S-matrix elements which constitute the expectation values ([ST)) . (|55|) . ([M)) . (|100p 
and (jlOip by using the cubic interaction (J116p . We here only consider the first order interaction and do not deal with 
loop diagram like FIGl3] Under these assumptions, we consider the S-matrix elements like 



{(p particle x 2; graviton x 1 out|0 in) . 

Then, we only consider ([57]) and (|55|) because the others do not include such S-matrix elements. 
Now, we derive the Feynman rules. Recall that 4>im and (f>io are written as 



(/.£„ = j di^-7^ {'J'tArn + -^tidtin) , 4>m 



duj 






(120) 



(121) 



where uj^g satisfy Eq. (j47p and the boundary condition in "in- region" (j73p . Then, we can deduce the Feynman rules 
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as follows: 



'l^e'm'lbtlm) = -T^^te^U'^rr^m' ' </*/„' M^fm) = -J^^te^U' ^rnm' > ^t^i' o\<te) = ~7^^'^i^U' ' (122) 



We can also obtain the Feynman rules for ^^ and ip^-^ . Note that i/jf„^ and V'm can be expanded as 

VL = / du;^ (wf.e™ + "f^rffL) . V'fo = / rf^4= (^f.cf, + u^^c^j) , (123) 



(-) \ LA/t LUCin LUC UJVtn I ' ' iU I /(-. 

where u^g satisfy the Zerilli equation and the boundary condition in "in-region" (j73p . Then, the Feynman rules for 
i)^ become 



The same rules apply to t^j^'^ . Using such rules and Eq. (|118p . we can perform the following calculations 



^^ntM5£™,6^^,,„M&^'^",n") =H\ ^'^* '^*^»"*'^tm2\f3m3 ^fiL>f2m2 <^f3m3 Mf £m:^l',„' :&^/',„" ) 



int 

= CS",,„ „ y I Ar*dt ^^^^^^ ul, (ut.ut,,„+4 -«t ») . (125) 

£ m ;£ m / ^ \ /o, „ ,', ," "'-V"^"^ I m t m J ^ ' 

Note that only one term in (|118p is related to this calculation. Consider another state Mf^mj^ /„/ ',c „ „). It is 
easy to obtain 



^^ 



H.nAd'.£„.bl,^,^,,c^^.,^„) = 5: / dr*dt/.„,Ct-^,^o<L^0.,™,0.3o|df,„„6^,,„,,i^,„„) 



int 

= C'r, .,^y fdr*dt^^^^u^Jut,,ut,,„+4 ,ut„ „) , (126) 

where the two terms in (|118l) are relevant. In the same manner, we can consider the terms related to ifj^'^ . For the 
state MX, e.'„,-^,y'."„/').-e have 



H.nAd5Z,bt,,,^„bt„,,,^„) - E/'^-*^*'^-<^"2^.3™3<™!^^-=^^3™3l«>l^„.6^,<<„,^; 



-^^ 



int 

= D^r,,, „ y f dr*dt- ^^^^^^ u^Y (ut.ut,y, - 4 ,4, ,) . (127) 

£ m ;£ m / ^ I /Z ~' ~ ^£ \ J £' ui £ £ m I m I ^ ' 

Note that additional sign comes from the antisymmetric property ^^^^'.(."^" = ^-^^p^^" -p'^' ■ These results suggest 

h. 



^£m 

m \£ m £ m \£ m 

that it is useful to define 



if,",v.,.„,„ = y^ I dr*dt^=^^^^==ut;[ur,,u%„„+u'':„,u 






TTUj£ _ 

ui'£'-u:"£" — 



y f dr*dt^^^^u^^* (u'^*,ut., -u'^lyut.,) . (128) 
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In fact, using these notations, all S-matrix elements can be written as (coefficient of <\118^ )xK, , ,, „* or (coefficient 
of (fTTOl) )xg"5., „„„*. Noticing that the effect of the first term of (fTUgi) on i7"/„, „ „„ is given by 



I 



dr*dt- 



f 



,RW* 



Using these notations, we calculate (|97| and ((98|) in the following 






U,.,U,, 



(129) 



1. Calculation of \98\ 



1 and one graviton " or 



First, we consider 7/'£„. In the formula (1981) . we only need the cases "n, ,„ = n, „„ 
"n , = 0, n ,, = 2 and one graviton". From (jllSp . for both cases, this graviton should be "d-particle" . In the 
former case, (|98|) becomes 



^EE E ^(0-1^^'..'^' 



L2 



, ,d ,„,,„ ,„ out 



j'" £'" 0<m"'<^"' 



^^'"EE^ 

i'" ni'" 



)(f6'^,, ) ,d^,„„,„ ,„ out 10 in) 






UJ £;uj £ UJ i-uj £ 



^£ m jy'uj £ 
r^£m;£m ^ uj' £;uj' t 



(130) 



where we have introduced the regularization (6 , „ ,b ,, „ ,d^, ,,„,,, ,,, out|6 , „ ,6 ,, „ ,d^,„„,,, ,,, out) = (Siio 

^ \ M iJTn ' /,! fym ' I.} i^ vn ^ UJ £m UJ tm UJ ii m ' 



u) £m UJ £ra uj £ ni 



0))^ = L^. Here, L has the dimension of the length. Note that we used Cf,™ , .„ „ ex (5„, ^' , ^" in the last equality. 

In the above calculation, we used the formula J2uj"' = L j du) . We can get the same result as this for the latter case 
and adding these results gives 



^ '' T j / J lin^lvn I UJ t\bj i UJ t\uj £ to t.uj i uj i\uj t\ 



(131) 



We only calculated the case that (^^^^ is "b-particle" so far. The same calculation leads to the result for "c-particle" 



]^ /* I III 1 r III III III III 

- \ da;"'2_]|Q''o/o Y^J £luj" £ ^Z' t,uj' £ 



UJ t.UJ i UJ 



' t.u" l\ 



(132) 



and that for "d-particle" 



T j / J tm,£m ^ i^ l-^uj £ w £:uj £ w l-.uj £ w £;w £j 



L 



rf-"'E|^. 



£ -2m 
m:i—m 



Ik'^'!''^'",, *K'^',"^"! +K'^„^ „ *K'^, " n 
Y^j£-,uj £ ^^uj £;uj£^ -^^ ■•'•■■• ' "./'>•.." 



UJ £\UJ £ UJ 



y: \ 

£;J £\ 



(133) 



Taking the summation of I and m in ([55]) . the effect of ([M]) on (^5]) is 



i/rf-"'E E E|<'^^r{C'Z'' 

■^ n n^ ,.,^ ^ n r.iii 



UJ i\UJ i UJ t\UJ i UJ 



£;uj £j 



i -e<m<f 

Next, we consider ip^'^ . Same calculation leads to almost the same result as (|134p 

^ f^,,'"\^ \^ Y^ I n^"'2"i (ttuj'"/" *ttJ"£" X ttJ" I" * TTUj'" e'" \ 

— / aUJ / / / \J--'fm-fm \ tl I I, II n n I. I . ~\- £1 II. II I, tl I. II. > 

•^ ^ -£<m<£ £'" 



(134) 



However, because of H'^ 



„ii ~ —H'^n„ii i„i, this vanishes. 

£ UJ i -.uj £ ^ 
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2. Calculation of {5 



First, we consider ih'^ . In this case, we must consider n, , „ — \,2. li n, , „ —2, the same calculation as before 
yields 






Z 

cj £ rn 



out) 



1 



= 7/'^-"'ER™; 



2m 






(135) 



Next, wc consider n , ~ 1. In this case, we must consider "another i/i-particle and one graviton". Furthermore, 

if this additional (/)-particle is "b-particle" with {ui,£,fh), we must consider a different mode {uj ,£,m) ^ {u),l,ifi). 
Then, we have to calculate the foUowings 

I 2 
^ I 2 

E ^"^ ™l^!'fm' 4,v V-^graviton x 1 out) , 



E ko in|6^^„^,4,-^, V'^'graviton x 1 out) 



(136) 



Note that the rule ^^ ~ L J dut should be used when we evaluate these quantities. Taking look at (|118p . we see that 
the graviton in the first equation of (|136p must be "d-particle" . The result is given by 

Lfdu;'" y y |(0in|6^, ,6^. ,d^,„, out) 



{uj,i,m)^{uj ,^,m) £ rn 

L I duj'" y y Icf "+" 

Z ^ Z ^ £m-,im 

{Q,e,m)^{uj' ,i,m) i'" 






(137) 



We can also calculate the second equation of p36p . In this case, the graviton must be "d-particle". The answer 
becomes 



L du:-J2J2\c£- 



K 



at,J e. 



(138) 



In the last equation of (|136p . we have to consider all types of gravitons. More precisely, the sign oim — m determines 
"particle" species. Thus, we obtain 



^/^-"'EEfE{ 



UJ, £,771 t l-m' 



(Oinlb"*,, ,d*.y,„, out) + (Oinlb"*,, ,d^. ,d^,„,,„ ,„ out) 

^ I /,) Pm^ lilPm.^ UJ P. 771 ' ^ I /,) Pm^ iilPm.^ LO P 771 ' 



UJ im'' uolm'' UJ p. 771 



(0in|6"^,, ,d't.,c^,„,,„ out) 



u),i,in l'" 



Eh 



if", " - 






if", " - 

u) 1-Sj£ 



^em;£-m 



K 



UJ t,a£ 



(139) 



This equation is divided into three parts, but only one has nonzero value for fixed m and m because C^V}/ , ^,, ,, oc 

^ r- ; J ell m :£ m 



m.m -\-m 



Then, after taking the summation of m. , we have 



Ljdj"yy\ci 

LU,P,ni £'" 



m—m 
fh:P77i 



2 









(140) 



Summing up Eqs. (|135p . (|137|) . p38p and (jl40p . we see that the contribution from the tree diagrams to ([W|) becomes 



dLj'"y y y c,^_7+^ 

ujI -l<fh<I P'" 



2 


III rrr 




K"'/, 




LJ i-.ul 



(141) 
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So far, we have only considered the case there are at least one "b-particlc" for <j>. Calculating "c-particle" and 
"d- particle" and taking the summation of £ and m in (|95p . we obtain 



duj"'y V V V c'-7+^ 

^ i.i'" I -i<m<i -I<fa<l 



2 /// /// 2 



(142) 



Next, we consider ij}^^ . With similar calculations with ([M]) . we can get the result by replacing C with D and K 
with H in p42p . The resultant expression becomes 



" £/" ,t-t<rii<l-f.<m<f. 



H'^, 



e-sii 



(143) 



3. Summary 

To summarize, substituting these S-matrix into (|95p . we can write down the effect of interaction on Hawking 
radiation from the tree diagrams as 

I /// 1 2 / /'/ '" '" /// '/' /'/ '" " 

^l 2m / T^u I * jj'tj £ I Tfui I * T^ui I 

X |Wm;£m| \-^ uj' l,uj" I ^ u' l-uj' I '^ ^ u" l,uj" I. ^ u' t.u' 
Q.u)' i'" I-i<fa<l I 



ini-£rrt 






t m+ra 
£m;£rn 


2 




111 



(144) 



where a and /3 are calculated in ([75]). Recall that we should use the rule ^^ = L J duj when we evaluate the second 
term in the above formula (|144p and ^^' in the second term should be j- J duj as understood from (|93|) . 



D. Cancelation of Divergences at horizon 

There are coupling function hint which is proportional to l//(r) in the definition of i^ /.. //.//, so it seems that 
K'^f^ ,, „,, would not converge. However, as you can see, the divergences cancel out each other. In the interaction 
(|113p . terms proportional to l//(r) are 



or 



A + 1 



^r-V^fm + ^dr'tpf,n <i>i' „^' <Pl" m" 



2t 



Isf 



2r , 



rf \ 



^dr4L - ^^f™ ) '^i'rn'dr'^r 
, 7s/ Isf I 



A + 1 



2r 



52 „i*Zl 



/ -^* J. 



7s/ 



.11 

Isf 



+ ^{ 9r-^f^^l'„,4t"„^" + dr'i^f^dr'(t)i' ^> dr'(J3i- ,n" " 2 V'f„*>£'™' ^r- </'/',„" 



(145) 
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Then, the contribution from these interactions to K'^i„, ,, „i, becomes 



dr*dt- 



\/& 



iWOJ UJ 



I o2 Z* ■ it>* ■ 0* o2 Z* ^ 4>* a <P* 



where w^ and u"^ satisfy Zerilh equation 



rf 

H — — 

Isf 



,(146) 



and equation of motion for (jj 






In these equation. Vzeriiu and V^ dump sufficiently first as r* — > — oo (or near horizon). Therefore, u^ and u'^ behave 
as 

ut,^a{t-r*) + i3{t + r*) ' ^'^'> 

Substituting these expressions into (|146p . we see that the contribution from near horizon vanishes. Indeed, the 
contribution from the pattern {u^,u'^,u'^) ~ (a, a, a) is clearly zero and that from the pattern (a, a,/3) is cancelled 
out by the pattern (/3,/3, a). Thus, it tmms out that K'^f , „,, converge by taking into account the equations of 

motion for -0'^ and (j). 

E. Deviation from Planck Spectrum 

Now, we estimate the effect of interaction (I144D . To begin with, we estimate K'^f'., ,, „,, and iJ"/., „,„ defined by 



T^uit \ ^ / J * ji '^int z* ( 0* -0* I -0* 0* \ 

K , ,, n 1,11 ^ y / dr dt , ^=u, ,« [u^,„,u^„„„ +u^ ,„,u^ „ „,, , 






UJ I ■,u I / y I /q ', ," "' \ u: I bJ I U) I UJ " 



where u^-^^'^ satisfy the boundary condition of "in region" (j73p . Using the geometric optics approximation, we can 
deduce the asymptotic solution 

"'^^ ~ \ e-*'^" - e-'^'^""" (as u ^ oo) ' ^ ' 

where u = t — r* and v = t + r* . Then, in order to estimate the effect of interaction, we approximate the solution in 
two ranges as 

z,w,0^re— -e— (for -c^<«<-iln«6) 

(for - i In Ke < u < oo) 



and we call the former range "outer area" and the latter range "inner area" (see FIGU]). 

We consider the case that all u^^^^''t> in K^f„, „ „,, and H'^f„, ,, „„ are e"'"". For both "inner area" and "outer 

UJ a -.u I UJ p \uj t 

area" , this case does not contribute to the integral because this pattern is proportional to the delta function 

^^^iuj+uj'^J'Y oc 5(w + w' + w") = . (150) 
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u = — ln{Ke)/K 



FIG. 4: We separately consider near u ~ oo and the other area. We put this boundary at u — —^ Iukc. 

Since we only consider the positive frequency mode, the delta function gives zero. Hence, we need to consider only 
the mixed cases. We can show that the contributions from "outer area" always vanish. We have already considered 
the case that all u-^^^^'"^ are e~*"'". Let us consider the cases that one of u is e~'"". If we consider only the "outer 
area" , this effect is negligible because the integral becomes 



- 6{u; + u;' +uj") / dr*/i„t(r)e*('^+"' """)'■* =0 . (151) 



Similarly, we can show that the contributions from "outer area" for other cases are also zero. Therefore, we can 
concentrate on the contribution from "inner area" . 

We consider only the dominant terms in interactions (|109p and (|113p . First, we check the interaction terms 
proportional to l//(r). However, these terms do not affect the integral. Because we assumed (|149p . the same 
calculation as the last subsection shows that these terms do not affect. Secondly, we consider the terms which are 
proportional to / x 1// in pi3p . These terms are 



— (dl,ll^^*4>dr-(t>- dr'1p^*<j)<j)^ - — (dr'll}^*(l)dr*(t>- i^^* H) , (152) 



again these terms do not affect the integral from the same consideration as the last subsection. 
There are 0(/°) in (fTT^ : 

^(r)9^.V'^*0(^ + C(r)a^.V'^*9r*?i9r*0 + ^(O^^>c'^*?^ + 5WV'^*00 + £'(^)V'^*9r*0a^*'^ • (153) 

And also exist in ([TM]); 

-xl;^^*^dr.c^ . (154) 

r 

In (|153p . we can show that the first three terms becomes zero by using the same method as the last subsection and 
(jllSp . Furthermore, the other terms in (|153p behave B(r) — > and D{r) — ^ as r — > oo; and coupling function in 
(|154p dumps as r — >■ oo. Then, these terms do not affect in the asymptotic region. This means that we can estimate 
the effect of these terms using not coupling functions but their value at the horizon , that is B{r = 2M) ^ 1/M, 
D ^ 1/M and 1/r ^ 1/M. Other terms can be neglected because their coupling functions are proportional to f/r. 
In fact, these terms do not affect both near the horizon and in the asymptotic region. Therefore, summarizing these, 
the dominant interactions in "inner area" are given by 



M V 



r*0(/) + V'5r- 05^*0^ (155) 
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and 



Thus, K'^f., , n ,11 can be estimated as 



K, , „ 



dr*dt- 



M 



1 1 



1{^«^*^9..</,} 



Note that we ignored the constant factor. And i/"/„/ ,/ „,/ is 

1 1 



^^ UJ i UJ £ Lot ' cj £ ' LJ £ 



^uj' £'-u" £" 



dr*dt- 



Vi 



UJU) U! 



' M 






0* 



(156) 



(157) 



(158) 



jZ.(p p — iujv 



However, from (11491) . we see this vanishes. Then, we consider K'^f„, „ „,, hereafter. We assume 

' ^ '^ ' u £ \uj £ 

e"'^*^ . In the above, the pattern {u^ ,u'^,u'^) = (e~*"'',e~*'^'', e"'"^") has been ahcady considered. Hence, 
at least, one of {u^,u'^,u'^) should be e~'^^^ . Among these cases, the most dominant one is {u^,u'^,u'^) = 
(e-'^<=""",e-^"'',e-'"''). This is because u e {-^lnKe,oo) implies that 



- 0(e) . 

Thus, the case with no derivative of e"'^*^ is dominant. From these consideration, K'- 
follows; 



(159) 
can be estimated as 



uj £ -u £ 



dr*dt- 



1 1 



^i{uj +ti) )v 



V ww'w" M 

UJ LU 1 



UJ M 



dve^^^ +" )" 



due^ 



— In K,e 



Uj' Lj" 1 1 



kU 



(160) 



The integral in the above equation diverges. However, this divergence occurs because we only considered outside of 
the horizon. In fact, since the horizon fluctuates due to metric perturbations, we must extend the range of integral 
to inside the horizon. We assume the range can be extended by e. Then, the integral converges as 



KU K 



So, K^f.i, II .„ can be estimated as 



l\ I nl II „ll 



yujbj'uj" 



(161) 



(162) 



Finally, using this result, we can evaluate the nonlinear effect. Assuming C^,™ , „ „ ^ 0(1), (|134p can be deduced 



as 



L J ^ K^AP Oj' +L0" ^ L K^M"^ w' + w" ' 

where LOcut is a cut-off frequency. In the same manner, (|142p can be estimated as 



7 h'-T. 



e^w (jj OJ 



IE- 



(163) 



(164) 



LJ""^ k'^M'^ iJ +ujY L^ k^AP (uj' + w)2 
Here, we note the Bogolubov coefficients are explicitly given by 

/3„,,„., i_^'-«-/«g-W{2K)r ( 1 + !^ ) , (165) 



1 



VKg7rc^/(2K)p ( 1 + _ 
K 
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from which wc obtain 



J2-KLJ / K, 



1/3. 



1/3. 



1 



1 



Now, we can calculate the contribution from (|142p to ((95)) as 



EE 



(eW« + i)|;3^^^^,^,|2^ 



M2L (w' + Cjf 



i-KKUj' e^'"'^/'^ - 1 



eV^,, e^-Z^ + l 
K3M2i e^'^'^/" - 1 



(166) 



(167) 



And. the contribution from (|134p becomes 



duj'duj'' 



K 






eV„. g^./K + g-TT./K g2^3 g2,r./K^j_ 



K^M'^L €''"/« - e-™/« K^M'^L e^'^"/" - 1 



(168) 



Taking into account that u^ has dimension of the length, we can recover £p. Thus, it turns out that the cubic 
interaction modifies Hawking radiation by the term proportional to 






coth(7ru;/K) 



(169) 



This effect is suppressed by the Planck scale. Therefore, the deviation from Planck distribution is very small but 
always exists . 

V. CONCLUSION 

In this paper, we have studied Hawking radiation from fluctuating black holes. For this purpose, we constructed 
the canonical quadratic action for metric perturbations and also obtained the cubic interaction terms in 4-dimensional 
Schwarzschild background. Taking into account cubic interaction terms in the action for a real scalar field, we derived 
the deviation from the Planck spectrum by considering only tree level diagrams. It turned out that the deviation is 
proportional to coth(7rw/K) . This result may have implications for black holes produced at the LHC and in the early 
universe. As is mentioned in the introduction, our analysis is also related to the trans-Planckian problem in Hawking 
radiation. 

There are many works to be done. First of all, it is intriguing to investigate higher order corrections for metric 
perturbations, i.e. self-interaction terms of ip- As the amplitude of quantum fluctuations of the metric is estimated as 
ip/fh, these higher order corrections must be important if black holes become as small as Planck scale. Considering 
these terms, we can investigate "non-Planckianity" of Hawking radiation from black holes like Non-Gaussianity 
in inflation [l^. Furthermore, from the AdS/CFT correspondence point of view, these interaction terms will be 
important for three point correlation functions in CFT. It is also interesting to generalize the present analysis to 
higher dimensional Schwarzschild black holes. This generalization is important because higher dimensional black 
holes might be created at the LHC if dimension of space-time is more than six. The calculations would be more 
complicated, but the linear equations of motion have been already derived in [1^. Therefore, it would be possible 
to compute the effect of interaction even for higher dimensional black holes. It is also possible to consider gauge 
flelds and fermions instead of a scalar fleld. Finally, since Einstein theory may not be valid near the Planck scale, it 
would be important to change the gravitational theory from Einstein theory to other theories motivated by quantum 
gravity. In fact, string theory predicts higher curvature corrections in addition to Einstein- Hilbert action jl7| . Again, 
the linear analysis has been already done in [iSl . Il9l | . Hence, we can extend our analysis to these cases. 
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Appendix A: Scalar Harmonics and Vector Harmonics 



We classify metric perturbations into scalar and vector perturbations in terms of the symmetry in 2-dimensional 
maximally symmetric space with the line element ds^ — "fabdx'^dx^ and the curvature k = ±1,0. They can be 
expanded by scalar harmonics and vector harmonics. In this Appendix, we summarize our conventions. 

Scalar harmonics Yfe are defined by the eigenvalue equations 



^ r. \n — Is ^ p- 



'p w 



(Al) 



where the eigenvalue 7s depends on the index p. For example, if the 2-dimensional curvature is positive k = 1, scalar 
harmonics are nothing but spherical harmonics Yfm and the eigenvalue is given by 7s — i{i + !)• In this case, p 
corresponds to £m. We impose normalization conditions 



jd^x^Y.Y^, 



Using scalar harmonics, we can define vector harmonics 



V, 



a ^ah 



^a^y 



p\h 



(A2) 



(A3) 



where tab = \/l£ab and £23 = ^1 = ^£32- The vector harmonics satisfy the transverse conditions V^^u = and 
eigenvalue equations 



Note that the vector harmonics satisfy the normalization conditions 



d^Xy/^VpaVp'''^^s5 



p,p 



(A4) 



(A5) 



Appendix B: Quadratic Action in (A)dS-Schwarzschild background 

1. 3+1 decomposition 

In this Appendix, we calculate the quadratic action for metric perturbations using 3+1 decomposition. The back- 
ground spacetime we consider is a (A)dS-Schwarzschild black hole 



ds^ 



-f{r)dt' + l/f{r)dr' + r'^abdx'^dx' , 



(Bl) 



where jab denotes the metric of 2-dimensional maximally symmetric space with the curvature k = ±1,0. Using this 
metric ansatz, Einstein equations with a cosmological constant A i?^i/ — Rg^^/2 + Ag^i, = read 



f =ik- f)/r - rA 



I 



-2(A:-/)/r2 



This can be easily solved as f{r) = k — 2M/r + Ar^/3. Here, M is the constant of integration. 
The lapse function N and shift vector N'^ of this space-time are given by 



^ - // , ^' = 



The induced metric is 



fill 



9i] 



V 



\ 



r'^lab 



Apparently, the extrinsic curvature of a t = const, surface vanishes 



(B2) 



(B3) 



(B4) 



K,, = 



(B5) 
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The spatial curvatures arc calculated as 



Rrr^-^, Rab^U-f-^Uab, R = 2A . (B6) 

According to 3+1 decomposition, Einstein-Hilbcrt action becomes 

S = dt(fxN^ (r + K'^K^j - K^ - 2A') , (B7) 

where g is the determinant of 3 dimensional metric and K is the trace part of extrinsic curvature. From now on, we 
call R-2A "potential term" and K'^K.j - K"^ "kinetic term". 

Let us define g = Ny/g and L = R + K"^^ Kij — K^ — 2A. In order to calculate the quadratic action, we must expand 
g and L; g ~ g^°' + g'-^' + g^"^^ + • • • and L ~ L'-^^ + U^' + L^^> + ■ ■ •, where (n) means the n-th order perturbation. 
Substituting these expressions into (jB7[) . we obtain 

S = f did3x(ffW + g(i) + 5^'' + • • ■){L'^'^ + L^'^ + L^^^ + ■■■) 

In the above action ([B8|, we notice i^"' = from dBS]) and ([B6|. We also have ^'^'i'") + g'-^'^L'^^'^ = because this 
term is proportional to background equations of motion. Furthermore, because the extrinsic curvature of background 
vanishes, L'^^ is given by the first order perturbation of R, M^\ and L'^^-' can be calculated as 

i(2) ^ ^(2) + ^(1) vj^^H)^^ _ (^K^y (B9) 

Hence, the action becomes 

S = j dtSx{N^)^^^ R^^^ + j dtd^xN^iR^"^^ +K^^^'^K^^\j-{K^^^)A+--- . (BIO) 

This tells us that we only need the first order lapse function, shift vector and extrinsic curvature. Hereafter, we omit 
" " ", that is, R\, means first order perturbation of 3-dimensional Ricci tensor, for example. 

Finally, we list up formulae for calculating perturbed quantities. The 3-dimensional background and perturbed 
metric are gij and hij , respectively. With this notation, the first and second order perturbations of the inverse metric 
are 

g(i)u ^ -g^'^g^'hki , 

g'^'' = -.9(')'V/iM-ff'"/"ff^'/i™«/ifei . (Bll) 

We can calculate Christoffel symbols from the following formulae 

r^% = Ig^'^'" {hir,k + hik;, h,,,i) = g^'^h^iT^^^ . (B12) 



And, Ricci tensor are given by 

7?(i).. _ r(i)' _r(i)' 

R^'\, = r%, - Ti'^U, + T^%r^'^Z - r(^H™r(i)^ . (bis) 

Note that ; means covariant derivative with respect to the 3-dimensional background metric g^j . 
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2. Scalar Type Perturbations 

We treat scalar perturbations using Zerilli gauge. Later, we will change the gauge to a more convenient gauge. In 
the Zerilli gauge, metric perturbations are expressed by 



(B14) 



( fH H, 
Hi H/f 


\ 






\ 


r^K-fab 
J 



Then, up to the first order, the lapse function and the shift vector are given by 

H 



^^ = V/ ( 1 - Y ) , N'-^fHi, TV^ - , Nr^Hi , Na^O 



and the induced metric becomes 



hij — 



f H/f 




V 







\ 

r^K-fab 

J 



(B15) 



(B16) 



First, we calculate the second order perturbation of "potential term" in the Einstein- Hilbert action. Using formulae 
(jBlip ^ (jB13|) . in the Zerilli gauge, perturbed quantities are calculated as follows: 



• inverse metric 



-,(!)"■ = _fH , £,(1)™ = , 5(1)°'' = -^AV'' 



(B17) 



^(2)r. ^ ^^2 ^ gi2)ra 



, 5(2)-^ = A.K^r' 



(B18) 



• Christoffcl symbols 

p(l)r 

rr 

p(l)a 






(B19) 



<2)r _ 



HH 



r^'Ka - -^ , r^'Kb = -fH (rH lir'K)'^ 7a. , 



r^'^". = tAiKH^^ , r(2)«, = -l-KK's^, , r^'^t, = -^ f a.a" + a>^c - k^^ 



2r2/ 



7c6 



• Ricci tensor 

^^ ra 



k"-^-k'-C-k' + --^M''w 



2/ 



2r2/ 



2r 2 



:H 



\ab 



f[rH--{r'K)' 



UrH-'^ir^K) 



-rfH - -A-l'',, 



(B20) 



lab . (B21) 
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Scalar curvature 






l"^ ' J.2 



„2 f' 



2kK + 2rHf -r^f K + 2fH + 2fr{H - 3K - rK 



(B22) 



NV9 



(VV^)W=.2V^(^A-+|-|) 



(B23) 



Using these results, we can obtain "potential terms" of the action . The second order scalar curvature perturbation 
can be calculated as R'^^^ ~ g'-^^'^^ Rij + g'^^^^^ R'^'^^j + g"^^ R^-' . The first two term can be easily calculated by using the 
previous results. Using Ea. (|B13|) . the last term can be written as 



9 



ij (r'-^)l _r(2)i , p(l)i p(l)m _p(l); p(l)r 

I ij-.l il-.j ^^ ij Im 1771 j 



(B24) 



In the above Eq. (jB24[ ). the last two terms must be calculated using the formulae (|B19p . while the first two terms can 
be calculated by using the integration by parts 



^N {g^^T(%, - ff^^r(2)^,^) = f M-diN) (g^=T^)\^ - 5^'r(2) 



V5(-5.A^)(^r(')a.-/r('^™ 



(B25) 



Thus, we can show that the second order "potential terms" become 



VI 



HH\''\a - r^/HK + -r^fK^ - r^fKH + - (H - H) r'^R^'^'> 



(B26) 



Next, we calculate the "kinetic terms" of the Einstein-Hilbert action. Here, we define Eij — NKij. At the first 
order, we have 



Er 



1 



H - 2fH, -f HA , E, 



1 



Hl\a , Eab 



1 



2/ V" '■'"' " "V ' "'" 2 

Then, the trace part of the extrinsic curvature is given by 

E = g''E,, = 1{h- 2fH[ - f'H, 



r^K 



2rfHi) Jab 



2f 
K -—Hi 
r 



(B27) 



(B28) 



and the quadratic part reads 



E,,E'^ ^-^[h-2JH[-!'Hi^ 



1 
2^ 



^H^aHi^^ + _i_ (r^X - 2rfHi 



Therefore, the second order perturbation of "kinetic terms" can be deduced as 



(B29) 



V7 



1 



,1", 



K 



2 2 

--HiHi^''la-^K^-jHk + 2rHHi+r^kH[ . 



r^fHi 



Then, we obtain the total quadratic action for perturbed metric 



V7[^^^'"|a - r^f'HK + \r\fK'^ - r^fKH + ^ (H - H) r^R'^^^ 

-iifii/il"|, - ^k' - jHK + 2rHHi + r'kH[ + j {r^fHi 



(B30) 



(B31) 



We can derive the constraint equations taking the variation of this second order action with respect to Hi and 
H. However, we can not easily solve resultant constraint equations. Therefore, we need to change the gauge from 
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Zcrilli gauge to the convenient gauge defined by ([7]) . This gauge transformation can be obtained from ([6]) . The gauge 
transformation we use is 



sym H/f 



sym sym 
\ sym sym 



\ 





y 






V 



\ 




r'K^lab 



(B32) 



/ 



where the index Z means Zerilli gauge. This gauge transformation can be obtained by setting ^^ ~ (0, —it;, 0,0). 
Then, from (j6]), the relation between these variables can be read off as 



H'' = H + f w , H^ = Hi-w , H^ ^ H - 2fw - f w , X^ = - 



2fw 



(B33) 



Substituting the transformation rules (jB33p into the action (|B3ip , we get the quadratic action for metric perturbations 
in the convenient gauge as 



+ lHlfH-2fw'-f'w-^w] \a-{2rfH) 



(^)'^2 kf^ I, 



-H' 



-WW' 



\a 



(B34) 



Now, we expand metric perturbations by scalar harmonics and only consider the real mode (for example, if fc = 1, 
TO = modes are real). Of course, other modes can be calculated and we can get almost the same result as real 
modes. After the expansion, the quadratic action becomes 



/ 



drdt 



T-IT2 



7« ,' 



-^Ht - iswHi + ^w' + 2rHHi - 7, 



2/ + r/ 
' 2r 



wH 



where 



or 



gi = 2rfH - 2jjw , q2 = Is { H + f w 



ifr)' 



2fw 



kf^ 



H'+ls^w'~^Hq2 



9i 



lsT[r) \ 2rf 



1 / ' , rf - 2/ 



T{r) = rf -2/ + 7, 



T{r) y^ ^' ' 2rf 

Therefore, the Hamiltonian constraint equation obtained by varying the action with respect to H gives 

92 = 
and the momentum constraint equation obtained by varying Hi is 

2rH 



H^ ^vo- 



ls 



(B35) 

(B36) 

(B37) 
(B38) 

(B39) 
(B40) 



Substituting (|B38|) . (|B39|) and (|B40|) into (|B35p . we get the quadratic action which is a functional of 91. Using a new 
variable <j) = qi/T{r) instead of 91, we obtain the quadratic action 



drdt 



^I^<*«'-5lll<*«'-^I^''-M*' 



(B41) 
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where we defined A = (7s — 2fc)/2 and Vz is given by 

/ 



Vz 



2 ('2 1 „3 f'Z 



8A(fc + A)^ + 4(fc + Xyrj + 2(2fc + A)rV + r' ! 



4:{r^f' -2rf + 2rX + 2rk) 

+4/2(2A + rf) - 4/ {2A(2fc + A) + 2{k + X)rf' + r'f'^} 

2/ A2(A + k)r^ + 2,\^Mr^ + (9Ar - ir^K)AP + 91/^ 
73" (rA + 3M)2 ■ 



Finally, defining 



^^^ 



A 



A + fc' 



and using the tortoise coordinate, we obtain the canonical quadratic action for scalar type perturbations 

dr*dt \{dti>^f - \{dr'4>^f - \yz(r){Vf 



(B42) 



(B43) 



(B44) 



3. Vector Type Perturbations 

Next, we derive the quadratic action for vector perturbations using Regge- Wheeler gauge 






Va 





Wa 


sym sym 





sym sym 


0, 



(B45) 



where Va and Wa satisfy transverse conditions Ua'" = 0,Wa'" = 0. In this gauge, the zeroth and first order lapse 
function and shift vector are given by 



7V=v/7, ^^=0, N'' = % 



and induced metric becomes 







hij = I sym 
sym 



Wa 








(B46) 



(B47) 



First, we consider the second order "potential term" in the Einstein-Hilbert action. In our gauge, perturbed 
quantities are calculated as follows; 



• inverse metric 



g{l)rr ^ Q ^ g(iy 



f...a 



W 



giDab ^ 



J2) 



'■'■ = ^-^W^Wa , g^^)™ = , (7(2)afc ^ / y^a^fe 



(B48) 
(B49) 



• Christoffel symbols 
F 



— U , i ,,Q — —W , i ah ~ ^ \Wa\b + Wb\a) , 



p(2)r 



r(l)l = 4r Iw"' + ^W'' 



J [ a' \ J a 

Wa\W + -—W 



2/ 



' '■'' ~ 2r2 V I'' b ) ^ ^ be- ^IbcW 



T^^^L^-^WjJwt^wi'] . F(2):, = — u;^7^,7,, 



(B50) 



r(2)-^ = , r(2)-, = Lw-ujt , f(2)«^ = --A^a (^^1^ + ^^1^ 



2r2' 



(B51) 
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• Ricci tensor 



f' f ' 



(B52) 



This proves that the first order scalar curvature i?^^-' vanishes. Therefore, in ()B10p . we only calculate / d^xNy/gR^^ 
for the "potential term" . This term is easily calculated using the technique (JB25P and the result is 



Next, we calculate the "kinetic term". The first order Eij = NK.ij are 



/ JipY' \J 1 / Jy* 



2 f «^a + ^Va - Waj , Eab = -- {va\b + Vb\a) 



These lead E = g^^ E,j = and 



E,jE'= = 



^ (^a + K - «:) (-" + y - v<^'^ + ^ («.„ + .,.) (.''I^ + ."l'^ 



Then, the second order "kinetic term" is given by 

^1 / . 2 '\ / • 2 



V^ 



va ^ >-'a i^a 

r 



- \ Wa + -Va -V„ I I W" + -W" - w" | - 



2/r2 



kVaV^+v'\''\aVb 



(B53) 



(B54) 



(B55) 



(B56) 



Note that we ignored total derivative terms. Thus, we have obtained the total quadratic action for vector type 
perturbations 



Vi 



I 2 ,\ . 2 

-j\Wa^ -Va ^ "a ) ( w" + -f " - u" 



2/r 



^ ^fcVaV'' + f'l'^lal'b) + ^2 (fe'^'^Wa + W^aW'^'V 



2r2 



(B57) 



Let us expand metric perturbations by vector harmonics V^ and only consider the real modes. The quadratic action 
becomes 



2(A + k) I drdt 



1 /. 2 A /. 2 

— w H — V — V ] \ w -\ — V — V 

2 V r I \ r 



A 



A/ 



2fr 



2 ^V 2 

1' — W 

2 ' 2r2 



(B58) 



where A = (7^ — 2fc)/2 and the overall factor comes from J y/^Vp*Vpa- Taking the variation of this action with respect 
to V, we obtain the constraint equation which cannot be solved easily. However, this difficulty can be circumvented 
by working in phase space. 
From the Lagrangian 



L = (2A + 2k) I dr 
we can define conjugate momentum p as 



rl 

„ , w H — V — V 
2 V r 



A 



A/ 



2 ^J 2 

V -w 

r 






XT 9 

p = -— = (2A + 2k){w + -V- v') 
ow r 



which leads to 



w 



P 
2X + 2k 



-V + V 
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The Hamiltonian is 

H ^ I drpw — L 



(2A + 2k) / dr 



1 



P 



2 V2A + 2fc 
In phase space, constraint equation can be written as 

2 P , f P 



A/^2 



2 P I f P 
r 2A + 2fe V2A + 2fc 



At; 
/r2 



2A + 2fc \2X + 2kJ f 



2Xv 



. 



It is easy to solve the above constraint as 



fr 



-P + P 



4A(A + fc) \r' 
Substituting this result into the Hamiltonian, we have reduced Hamiltonian 



H ={2X + 2k) / dr 



iCfri /\ / P 
2 2X X] \2X + 2k 



fr 



P 



A/ 



4A V2A + 2fc^ 
which is a functional of w and p. Let us make the following simple canonical transformation 

Q = 'P, P = w, H{w,p) = K{Q, P) , 
where K(Q,P) is a new Hamiltonian. Then, we have 
K{Q,P) = H{P,-Q) 





= (2A + 2fc) f dr 


Yi 


ifr)' f\ 
2X XJ 


{ Q \\ f^' 


( "^ ) 


1 

2^^fp2 

r^ 




\2X + 2kJ 4A 


\2X + 2k) 


, we define 















5K _ 4A(A + k)f 
JP ~ ^ 



P . 



Thus, the action becomes 



/ drdtPQ - I 



drdtPQ - / dtK{Q, P) = 



1 



2X + 2k 



drdt 



1 ifr) , / 



Finally, changing normalization by 



V' 



RW _ 



4A/^ I 2 2A ■ A 



--Q 



f^2 



^\Q'-^Q' 



2yjX{X + k) 
and using the tortoise coordinate, we obtain the canonical quadratic action for vector perturbation 



dr*dt 



where we have defined 



\ {d.r'^f - \ {dr^^'^'^f ivwM^^>^^ 



V'w = ^(-r-/ +2/ + 2A) . 



(B59) 



(B60) 



(B61) 



(B62) 



(B63) 



(B64) 



(B65) 



(B66) 



(B67) 



(B68) 
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Appendix C: Symmetry in Quadratic Action 

111 this Appendix, we identify the symmetry in the quadratic action. As shown in [20| , we can rewrite Vz and Vrw 






where 



Now we consider the quadratic action (jB44p and (jB67p . Using VF(r), this action can be rewritten as 



- / dr*dt 



ij^HiP 



(CI) 



(C2) 



(C3) 



where ■0 represents (-0 ,-0^*^) . Here, the operator H is defined by 



H 



dW 



dt + dt.-[ui + W'+a,— ] , 



(C4) 



where (T3 is a Pauli matrix. Then, if there exist operators which commutate with i7, the canonical action (jCSP has 
symmetry. In fact, there are two such operators defined by 

Ql = (Jl-dr' + (T2W , Q2 ^ CF2-dr' - (JlW , 

I I 

where ui and 1J2 are also Pauli matrices. These operators satisfy anti-commutation relation {Qi,(52}=0. Then we 
can conclude that the quadratic action in the (A)dS-Schwarzschild background has "N=2 supersymmetry" . 
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